Observing the decay of orbital angular momentum entanglement, 
through experimentally simulated turbulence 
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We study the evolution of an orbital angular momentum (OAM) entangled bipartite photonic 
state for the case where one of the photons propagates through Kolmogorov turbulence, using the 
concurrence as a measure of entanglement. Quantum state tomography was performed to reconstruct 
the two qubit density matrices for a range of scintillation strengths. Our results give the first direct 
experimental confirmation of the existing theories for decay of entanglement due to atmospheric 
turbulence. We also show how the modal scattering increases with increasing scintillation and we 
discuss the impact of the scale at which entanglement dissipates due to atmospheric turbulence on 
free-space quantum communication. 
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Laguerre- Gaussian (LG) modes have recently attracted 
much attention within the quantum information com- 
munity, mainly because their infinite-dimensional Hilbert 
space allows information processing with higher dimen- 
sional quantum states (qudits) [l| for use in higher di- 
mensional quantum key distribution Q and long-range 
quantum communication [3|, among others. An LG beam 
with azimuthal index £ carries an orbital angular momen- 
tum (OAM) of £h per photon 0, 

Quantum entanglement is an important resource for 
quantum information processing and quantum commu- 
nication, but suffers decay when encountering a noisy 
channel 0, 0, such as atmospheric turbulence. Various 
aspects of OAM modes propagating through turbulence 
have been considered theoretically , including the detec- 
tion probability of OAM modes [3, 8, 9], attenuation and 
crosstalk among multiple OAM channels [1 
cay of entanglement for bipartite qubits [11 
most of these studies [8|,l9|, UJj are based on a single phase 
screen approximation [3j, which is only valid in weak 
scintillation, the case of arbitrary scintillation strength, 
which requires a multiple phase screen approach, has also 
been considered with the aid of an infinitesimal propa- 
gation equation (IPE) [l2j. These investigations show 
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that OAM entanglement is more robust when the beam 
waist radius wo is small compared to the turbulence co- 
herence length, which is given by the Fried parameter [l3| 
r = 0.185(A 2 /C 2 z) 3/5 , where C 2 is the refractive index 
structure constant, z is the propagation distance and A 
is the wavelength. The studies also indicate that entan- 
glement decays comparatively slower for larger OAM. 

From an experimental point of view, the problem of at- 
mospherically induced decay of OAM entanglement has 
thus far received little attention. It was shown, using 
coincidence counts, that the number of entangled modes 
(the Shannon dimensionality) decreases with increasing 



scintillation [14J . In other studies [15|, LL6[ the crosstalk 
among OAM modes has been measured, using a single 
phase screen to simulate the turbulence experimentally. 
None of these experiments directly addressed the dissi- 
pation of entanglement due to atmospheric scintillation. 

In this letter we present the first experimental results 
on the decay of OAM entanglement between two qubits 
due to atmospheric turbulence. To allow quantitative 
comparison between our results and the existing theories 
that are based on [3], we consider only qubits and im- 
plement the turbulence as a phase-only distortion on a 
single thin phase screen, using a spatial light modulator 
(SLM). We compute the concurrence [17] directly from 
the density matrices obtained by employing full state to- 
mography [lij to observe the entanglement dissipation 
as a function of scintillation strength. 
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FIG. 1: Experimental setup used to detect the OAM eigen- 
states after SPDC (a). The plane of the crystal was imaged 
onto two separate SLMs where the modes were selected, as 
shown in (b) and (c). The SLM planes were then imaged to 
the inputs of single-mode fibers, connected to APDs. Ran- 
dom phase fluctuations were added to the mode on one of the 
SLMs to simulate Kolmogorov turbulence (c). 
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A diagram of our experimental setup is shown in 
Fig. QJa). A mode- locked laser source with a wavelength 
of 355 nm and an average power of 350 mW was used to 
pump a 0.5 mm-thick type I BBO crystal to produce non- 
collinear, degenerate photon pairs via spontaneous para- 
metric down- conversion (SPDC). Down- converted pho- 
ton pairs are naturally entangled in terms of the OAM 
basis due to the conservation of OAM in SPDC [19]. The 
plane of the crystal was imaged onto two separate SLMs 
in the signal and idler beams, respectively. The SLMs 
served to perform pro jective measurements for quantum 
state tomography [18| by selecting particular pairs of 
modes for detection [see Fig. [TJb)]. The atmospheric 
turbulence was simulated by adding random phase fluc- 
tuations to the phase function of one of the SLMs [see 
Fig. Hie)]. The SLM planes were re-imaged and coupled 
into single- mode fibers, which extract the pure Gaussian 
mode from the incident field. Avalanche photo diodes 
(APDs) that were connected to the fibers registered the 
photon pairs via a coincidence counter (CC). All mea- 
sured coincidence counts were accumulated over a 10 sec- 
ond integration time with a gating time of 10 ns. 

The random phase fluctuation 0(x,y) that was added 
to the one SLM was obtained by the method of filtering 
white Gaussian noise 2o|, [H ] : 



In normalized coordinates the LG modes are 
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Here ^ r_1 {-} is the two-dimensional inverse Fourier trans- 
form, k_L is the two-dimensional wave vector in the trans- 
verse Fourier domain, is the sampling interval in the 
frequency domain and x(k±) is a frequency domain zero- 
mean Gaussian pseudo-random complex function, obey- 
ing x*(k±) = x(— because 0(x,y) is real-valued. 
The phase power spectral density is related to the re- 
fractive index power spectral density through ^g(k±) = 
27r/co^ < ^ ) n(k±, 0), where ko is the wavenumber (2tt/X). 
For comparison with existing theories, we used the Kol- 
mogorov spectrum (k) = 0.033 C^" 11 / 3 (HQ and 
to ensure that the random phase functions can repro- 
duce the Kolmogorov structure function we added sub- 
grid sampling points to the Fourier domain representa- 
tion of the phase function [24]. Our experimentally sim- 
ulated scintillation strengths ranged from wq/vq = to 
4, in 21 increments. For each scintillation strength we 
performed 26 realizations and a state tomography [l8| 
was performed for each realization to reconstruct the bi- 
partite qubit density matrix. 

The entanglement is quantified by the concurrence [ItJ 
C(p) = max{0, y/Xi — ^/X^— \A3 — a/AI}, where the A n 's 
are the eigenvalues, in decreasing order, of the matrix 
p((J y (J y )p*((Jy <g) <J y ) . Here * represents the complex 
conjugate and a y is the Pauli ^/-matrix. 
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where Lp'(-) represents the generalized Laguerre polyno- 
mials with azimuthal index I and radial index p, <f> is the 
azimuthal angle, r = (x 2 + y 2 ) 1 ^ 2 /wo and t = z/zr, with 
zr being the Rayleigh range (ttWq/X). The normalization 
and Gouy phase factors are contained in 
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The full spatial entanglement of down-converted pho- 
tons is experimentally accessible, if all LG modes with 
arbitrary i and p indices are considered [25[. However, 
here only LG modes with p = are measured. This al- 
lows us to express the state of a down- converted photon 
pair as = X^ C ^K)I — where \£) represents the LG 
modes withp = and |q| 2 denotes the probability ampli- 
tude for detecting the signal and the idler photons in the 
state \£) and \ — t), respectively. Here we only consider 
qubits defined in terms of the OAM states with £ = ±1. 
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FIG. 2: Experimentally measured intensity distributions of 
LG beams for i — 0, 1, 3 and p = 0, via back-projection from 
the fiber onto the SLM (see text). First column without tur- 
bulence (u>o/ro = 0). Second column for wo/ro — 1 and the 
last column for wo/ro — 3. 

Any classical LG beam propagating through turbu- 
lence is distorted by random phase modulations that, af- 
ter some propagation, turn into intensity fluctuations. 
Examples of such intensity distortions are shown in 
Fig. [21 These images were obtained in the plane of the 
BBO crystal with back-projected illumination [I?} of the 
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SLM that contains the phase fluctuations in the experi- 
mental setup shown in Fig. Ufa). The distortions of the 
LG mode implies that the energy of the original LG mode 
is scattered into other LG modes. 

A similar scattering process occurs in the case of an en- 
tangled quantum state propagating through turbulence. 
Given an initial pure state after turbulence it 

becomes Y,mnP™>n\t>m){t>n\ (ignoring the p- 

index), where p mn are the density matrix elements. 
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FIG. 3: Mode scattering under the effect of turbulence given 
by the coincidence counts for simultaneous measurements of 
modes with azimuthal index £a in the signal beam and in 
the idler beam, (a) With no turbulence, only anti-correlated 
coincidences are observed. As the scintillation strength in- 
creases, the mode scattering becomes more pronounced: (b) 
wo/ro = 1, (c) wo/ro = 2 and (d) w /ro = 3. 

We have experimentally measured this modal scatter- 
ing for different scintillation strengths shown as modal 
diagrams in Fig. [3] Without scintillation the modal dia- 
gram displays a perfect correlation between £a and — £b, 
as seen from the strong diagonal in Fig. 02a), which is in- 
dicative of the entanglement of the down-converted pho- 
ton state. As the scintillation is increased, one can see 
from Figs. [3fb) to (d) that the components of the quan- 
tum state are gradually scattered away from the diago- 
nal, first to the neighboring locations and then to loca- 
tions further away. This indicates that the mode coupling 
is strongest for neighboring modes and decreases as the 
mode separation increases. 

Our main result is the experimental observation of 
the decay of entanglement as a function of the scintil- 
lation strength, given by wq/vq. In Fig. |4] we compare 
the experimentally obtained concurrence values with the 
theoretical concurrence curves computed within the sin- 
gle phase screen approximation, as done by Smith and 
Raymer (S&R) and with the IPE (Tg. Since the 
experiment imposes scintillation on only one of the two 



photons, the calculation of the two theoretical curves was 
also done for only one photon passing through turbu- 
lence. The experimental results were obtained from nu- 
merous quantum state tomography measurements, from 
which more than 500 individual density matrices were re- 
constructed (26 realizations for 21 points). These density 
matrices were used to compute the concurrence. Care 
was taken to remove any negative eigenvalues that occur 
for these density matrices |28j . Such negative eigenvalues 
are caused by measurement errors — fluctuations in the 
coincidence counts caused by fluctuations in the photon 
number statistics. Each of the points in Fig. |4] repre- 
sents an average of 26 such concurrence values and the 
error bars indicate the associated standard deviations. 
Note that the experimental values indicate a maximum 
concurrence without turbulence of about C = 0.8. As a 
result we adjusted the theoretical curves by an overall 
factor to provide the best fit with the experimental data. 
The lower concurrence is due to experimental imperfec- 
tions and photon statistics. 
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FIG. 4: Concurrence as a function of scintillation strength 
(wo/ro). The (blue) diamonds represent the average concur- 
rence obtained from the experimental measurements, together 
with their standard deviations, shown as error bars. The solid 
(red) line represents the theoretical curve obtained with the 
single phase screen approximation (weak scintillation) and the 
dashed (green) line is the theoretical curve obtained from the 
IPE. The inset shows the same graph on a logarithmic hori- 
zontal axis to emphasize the scale of entanglement decay. 



There is a good agreement between the experimental 
results and the S&R curve. This is the first experimen- 
tal confirmation of these theoretical results as far as we 
are aware. The experimental results deviate from the 
IPE curve, especially for larger values of wq/vq. This is 
not unexpected, in view of the fact that the single phase 
screen approximation is only valid in weak scintillation 
conditions. It should be noted that both the experimen- 
tal setup and the S&R curve represent a single phase 



screen approach, while the IPE does not make this ap- 
proximation. The decrease in entanglement (Fig. [4]) is 
consistent with the modal scattering observed in Fig. [3l 
The most pertinent point of the results shown in Fig. |4] 
is the fact that it provides experimental confirmation 
of the scale of entanglement decay. Viewing the curves 
on a logarithmic horizontal scale, as shown in the inset 
in Fig. 3J one finds that the decay of entanglement all 
happens within an order of magnitude around the point 
where wq/vq = 1. The fact that both theoretical curves 
and the experimental data share the same scale indicates 
that, whether or not one uses a single phase screen ap- 
proximation does not make any difference to that scale. 
Substituting the expression for the Fried parameter into 
wo/ro = 1, we derive the distance scale at which atmo- 
spherically induced decay of entanglement occurs: 
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For practical free-space quantum communication systems 
the distances between repeaters would need to be shorter 
than Ldec- Although this distance would often be shorter 
than the Rayleigh range [14[, in weak turbulence this 
distance could be larger than the Rayleigh range. 

Interestingly, the scale where wq/vq = 1 roughly coin- 
cides with the transition region between weak and strong 
scintillation. Scintillation is classified as either weak or 
strong based on the Rytov variance, which is defined as 
o\ = 1.23C^kl^ 6 z n ^ 6 . For plane waves, strong scintilla- 
tion is said to exist when a\ > 1 [23| and for Gaussian 
beams it exists when o\ > (t + l/tf^ 0. 

The diagram in Fig. [5] shows the different regions in 
terms of the Rytov variance as a function of the normal- 
ized propagation distance t. For any particular optical 
beam propagating through turbulence, the value of o\ 
is proportional to t 11 / 6 . Three such lines are shown in 
Fig. [5] for different turbulence strengths. One can see 
that these lines start off in the region of weak scintilla- 
tion at the bottom of the diagram and then move up to- 
ward the region of strong scintillation as the beams prop- 
agate further. Eventually these lines cross the boundary 
into the region of strong scintillation [cr\ = (t + 1/t) 5 / 6 ]. 
However, at the same time (for weak turbulence) or even 
before it (for strong turbulence) it also crosses the line 
where wq/vq = 1 (dashed line in Fig. [5]), which is the 
scale where entanglement decays. The dashed line is 
obtained by expressing the Rytov variance in terms of 
w /r , which leads to a\ = 1.637 t 5 ^(w /r ) 5/s . For 
y;o/ro = 1 we then find that a 2 R is proportional to t 5 ^ 6 . 

Based on the above argument it seems that one never 
reaches the strong scintillation region with a nonzero con- 
currence, regardless of the approximations used, which 
implies that the single phase screen approximation can 
be used for all situations. We must, however, emphasize 
that this conclusion is based on our current results, which 



only consider the case where only one photon propagates 
through turbulence and where the basis is restricted to 
\£\ = 1. When both photons propagate through turbu- 
lence or when \£\ ^> 1 the difference between the S&R 
and IPE curves may become bigger, which may cause the 
scales for IPE and S&R to deviate significantly. For such 
a case situations may exist where entanglement survives 
deeper into the strong scintillation region. 
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FIG. 5: A diagram showing scintillation strength in terms 
of a R as a function of t on a log-log scale. Weak scin- 
tillation lies toward the bottom and strong scintillation to- 
ward to top, with the boundary between these two regions 
shown for both plane waves (a R = 1) and Gaussian beams 
[°R = (t + l/£) 5/6 ]. The dashed line (w /r = 1) is roughly 
where the concurrence approaches zero and the three slanted 
colored lines indicate the increase in <j r with propagation 
distance for different strengths of turbulence (left to right: 
C^{10- 12 ,10- 14 ,10- 16 }m 2 / 3 ). 



In summary, we present the first experimental confir- 
mation of the theoretical decay of entanglement in at- 
mospheric turbulence, where only one photon propagates 
through turbulence and the qubit basis is restricted to an 
OAM of \£\ = 1. More work remains to be done: the be- 
havior for quantum states with larger OAM and the effect 
of different scenarios, such as when both photons propa- 
gate through turbulence, need further investigation. 

This work was done with funding support from the 
CSIR. 
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